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A. OPTICAL DESIGN 


I INTRODUCTION 

During the preliminary design stage of this Investigation, 
approximately twenty different types of lens configurations were 
selected for the initial study effort* These lens systems were 
typical representations of three basic categories of physical 
arrangements: 

]* Catadioptric - Consisting of both refracting and reflect- 
ing elements 

2. Reflecting - Consisting mainly of first surface reflect- 
ing elements 

3- Refracting - Consisting totally of refractive lens elements 

The method used for selecting the optimum solution contained five 
essentia] criteria: 

1* Image Quality - The blur circle diameters must not exceed 
0*150 mm across the field of view considering 80 % of the 
ava i 1 able energy* 

2. Obscuration - The lens system should exhibit a minimum 
amount of energy losses due to physical obstructions within 
the system* 

3* Distortion - The target objective for the system should be 

3 arc minutes, or less, anywhere in the image format. 

4. Collecting Area - The lens system should provide at least 
2 

18*3 cm of collecting area for on-axis targets. 
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5. Physical Restrictions: 

a) System Envelope - The maximum dimensions of the 
package should not exceed a cylindrical volume 
of six inches diameter by six inches in length. 

b) Weight - The design goal was given as five pounds 
or less. 

The following criteria were used as additional weighting para- 
meters in selecting a lens configuration which would possess desir- 
able attributes in terms of producing a suitable prototype model of 
the theoretical solution: 

1. Material Considerations - High ultra-violet transmitting 
glasses, of low cost, ready availability, and good working 
characteristics would be needed. 

2. System Sensitivity Factors - The final solution should poss- 
ess relatively low insensitivity with respect to nominal op- 
tical fabrication tolerances and individual mechanical seat- 
ing errors in the lens mount. 

3- Ghost Image Locations - A suitable design would possess a min- 
imum number of ghost images close to the final focal plane. 

As was reported on in the first preliminary report there were four- 
teen individual lens designs that were brought to a point of semi-final 
completion. From these, only one solution exhibited outstanding perform- 
ance characteristics in terms of the first five rating criteria. This 
particular system was a six element refracting design, operating at a 
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relative aperture of f/1.0 and possessing good image quality across 
the entire image format. The limiting factor in this design select- 
ion was the necessity of using rather high index materials. Since 
high refractive indexes usually indicates dense flint glasses of 
comparatively low ultra-violet transmission, or heavy lanthanum type 
crown materials which are expensive, it was realized that this partic- 
ular solution would require careful study to insure the best material 
choices. 

FINAL DESIGN PROCEDURES 

There were many factors governing the final selection of appropr- 
iate materials for this solution. Since the specifications require 
high transmission in the low ultra-violet spectral region, only those 
glass types In the Schott catalog which gave values of 30 % or greater, at 

400 NM were considered useable. Therefore, the initial effort consisted 

* 

of adjusting the nominal theoretical solution to accommodate the exact 
catalog values. Schott lists eighteen separate lanthanum type glasses. 

Of these only six possess refractive indexes of 1.7, or greater: LaK-8 
( 713556 ). LaK-10 ( 720504 ), LaKN-16 (73^517), LaK-17 (788505), LaKN-l8 
(729542) and LaKN-19 (755531)- The last two types are special order 
items only and should be avoided for volume fabrication systems. LaK-17 
contains thorium and is four times more expensive than more conventional 
lanthanum glasses. Therefore, it was eliminated. LaKN-l6 is twice as 
expensive and not manufactured as frequently as other glasses in this 
classification. Therefore, the final selection of LaK-8 was based on 
its lower cost, desirable Abbe dispersion value of 53-83 and its slight- 
ly better transmission value, over LaK-10, in the U-V region. Identical 
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procedures were followed in the selection of SF-1 and SF-19 for 
suitable materials. There is a wider choice of heavy flint mat~ 
erials, however, so these two glass types could be chosen quite 
selectively in terms of how closely they matched the theoretical 
glasses of the preliminary design solution. These flint materials, 
of course, were less expensive than the lanthanum glasses by a factor 
of about four times. However, there was no suitable low cost substi- 
tute for the lanthanum type crowns and they had to be used. 

Upon completion of suitable optical material selection, the pre- 
liminary theoretical design values were then bounded so as to squeeze 
the solution towards the actual catalog values of the selected glass 
types. The optimization program uses all the lens variables, such as: 
lens radii of curvature, element center thicknesses and air spaces. 

Both damped least squares and ortho normal matrix inversion techniques 
are utilized in this optimization program to achieve rap-id convergence 
of the final design solution. 

DETAILED TECHNICAL DISCUSSION 

As the optimization cycles progressed it became quite obvious that 
the image quality below 4A0 NM degraded very rapidly. This is quite 
typical when correction over a rather broad spectral range is attempted, 
particularly at fast relative apertures. Therefore, as mentioned in 
the earlier report, this obstacle was anticipated. There are very few 
optical materials which permit a lens system to be corrected down near 
^00 NM, especially when good quality is required in the red portion of 
the spectrum. The exact transition point could not be predicted prec- 

-if- 
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isely since the tolerable blur circle diameter of I 50 microtis 
is quite modest and it was hoped that even at k20 to 430 NM, a 
satisfactory performance level could be attained. The dominant influ- 
ence affecting the cut-off point for the low end of the spectrum was 
the fast operational speed of f/1.0. At this relative aperture, the 
lens radii are quite short which, in turn, produces higher order abei — 
ration contributions. These surface contributions must then be bal- 
anced by equally strong contributions of opposite sign, which produces 
strong zonal characteristics. 

The system was also severely restrained by the requirement on the 
long back focus value of 1.00 inch. Ordinarily, in such a fast system, 
this working distance would be of the order of 10 to 20Z maximum of 
the equivalent focal length, rather than the required 50^. As the de- 
sign progressed then, there was a continual struggle to improve the 
image quality across the field and yet maintain the desired 1.00 inch 
distance. 

A careful check was also maintained on the individual positive ele- 
ment edge thicknesses and negative element center thicknesses. Both 
these values had to be maximized to insure adequate ease in fabrication 
and sufficient strength in mechanical mounting to assure long life under 
operational conditions. The final lens prescription that evolved from 
this analysis is detailed in the Table I below; 
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TABLE 1 



LOW SCATTER LENS 

SYSTEM: 

CONSTRUCTION DETAILS 


SURFACE 

NUMBER 

RADII OF 
CURVATURE 

VERTEX GLASS 

SEPARATIONS MATERIALS 

CLEAR 
APERTURE 
RADI 1 

1 

Entrance Pup i 1 

-2.7000 


0.971 

2 

3.599 

0.540 

LaK-8 

1.147 

3 

-6.109 

0.004 


1 .092 

h 

1.844 

0.633 

LaK-8 

0.895 

5 

4.251 

0.156 


0.640 

6 

-2.061 

0.250 

SF -19 

0.636 

7 

0.999 

0.179 


0.553 

8 

Glare Stop 

0.070 


0.555 

9 

-2.680 

0.200 

SF-1 

0.557 

10 

1.995 

0.002 

- 

o.6;27 

1 i 

1.811 

0.640 

LaK-8 

0.643 

12 

-1 .236 

0.004 


0.695 

13 

1.509 

0.380 

LaK-8 

0.624 

14 

21.849 

0 

0 


0.557 

15 

Image Surface: 
Flat 




NOTE: 

In order to achieve high volume 

productive capacity. 

the lens 


radii were constrained 

to equal 

1,000 inch, or greater. This 


factor did not exert an unduly severe restraint upon the 
optical design. 
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The performance characteristics for this final optimized 

solution are shown in Fig. 1, 2 and 3* These are conventional 

ray intercept versus ray height in the entrance pupil plots, for 

each of three field angles.- Fig. 1. is the information for the 

on-axis analysis. Fig. 2. for half-way out in the field, and Fig. 

3. shows the image quality at the extreme edge of the field. In 

each set of these plots, there are three separate representations. 

The first, located at the top of the page and labeled SAG FAN Y, 

shows the Ay, or comat ic type errors, for the sagittal fan of energy. 

* 

The second plot, in the center of the page and labeled SAG FAN Z, 
gives the AZ, or skew errors, for the sagittal fan. These two plots 
are given only for the center to the edge of the pupil since the 
negative side of the data will be a mirror image of the positive side 
of the pupil aperture. The bottom plot provides the image errors for 
the meridional or tangential, energy fan. The total aperture, from a 
relative height of +1.0 to a relative height of -I.O, must be shown here 
since there is no inherent symmetry between the two halves of the pupil 
aperture. The main vertical divisions are the image errors directly, 
for any particular ray height in the entrance pupil, and are given in 
.010 mm steps. Three separate wavelengths are distinguished on each 
plot, indicated by L (for low index), M (for mean index) and H (for 
high index). In this instance, as mentioned previously, L = 0.55 mic- 
rons, M - 0.4935 microns, and H = 0.45 microns. 

An examination of Fig. 1. shows that essentially all the available 
energy Is within a blur circle of about 0.080 mm diameter. Similarly, 
Fig. 2. which is for a position half-way out in the field, shows all 

- 7 - 


SPACE OPTICS RESEARCH LABS. 



the energy is concentrated within a O.O85 mm diameter blur 
circle. Finally, Fig, 3., for the edge of the field, shows 
about the same degree of correction. 

All these curves show the effect of very strong zonal errors 
which is quite typical of very fast optical systems. This necessity 
for balancing very strong opposing aberrations is usually the limit- 
ing factor in the level of image quality attained. For this case, 
the burden upon the optical design was intensified due to the enhanc- 
ed working distance requirement and the extended spectral range cover- 
age into the blue region. 

FINAL PERFORMANCE ANALYSES 

The performance characteristics of this lens system in a more 
pictorial form in Plates #1 through # LVIl, These plots give the 

J 

exact physical shape of the blur circle pattern for each wavelength 
and each designated field angle. Seven discrete steps in the spectral 
range extending from O.AOO microns to O.6OO microns are given at six 
different field angles from 0.0 degrees out to 7. 50 degrees, half- 
angle. The representations for the low blue region of O.kO microns are 
included even though the lens system could not be corrected this far 
down into the blue region. 

All the plots are initially scaled for one inch equals 0.005 mm. 

However, if the energy plots exceed this value then the program auto- 
matically scales the data to keep the total plot on the page. The 
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new scale is printed out at the top left hand position on the 
sheet. Plate I, for example, shows a scale factor of one inch 
equals 0.03^630 mm. This is for a wavelength of 0.400 microns, 
at a field angle of 0.0 degrees, as shown. In this plot 100^ of the 
available energy will be within 0.303 mm. Of course, the blur 
circle diameter becomes substantially smaller, as a lower energy 
concentration is considered. For example, if 751 of the available 
energy is considered, then the blur circle diameter becomes about 
0. 140 mm. 


A total summary of all this data is given in Table II below. 

TABLE 1 1 . 

SPOT DIAGRAM ANALYSIS 

BLUR CIRCLE DIAMETER, mm. ]00^ ENERGY 

FIELD 


ANGLE;SINE X 

0.4 

0.45 

0.475 

0.500 

0.525 

0.550 

0.600 

0,000° 0.000 

0.303 

0.126 

0.088 

0,066 

0.058 

0.062 

0.063 

0.026105 

0.313 

0.133 

0.095 

0.089 

0.074 

0.074 

0.066 

0.052210 

0.332 

0.150 

0.111 

0.086 

0.072 

0.071 

0.091 

0,078316 

0.361 

0.146 

0.105 

0.109 

0.109 

0.107 

0.086 

0.101|^21 

0.351 

0.153 

0.132 

0.112 

0.099 

0.093 

0.101 

0.130526 

0.0348 

0.217 

0.175 

0.149 

0.138 

0.105 

0.106 


The above results show that the system is well corrected over 
the entire spectral range extending from 0.A5 to 0.60 microns, even 
when considering the full 100^ of the available energy. The perform- 
ance at O.llO microns is also within the allowed 0.150 mm. if about 
80 to sot of the available energy is considered. Therefore, in terms 
of the minimum weight package, the best available solution of suitable 
glass types, and considering reasonable restraints for manufacturing 
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ease, the final solution represents an optimum compromise of all 
the factors involved. 

During the actual fabrication procedure for this lens system 
tiiere would be a so-called melt adaptation required. This process 
IS necessary because the actual refractive index of the manufactured 
glass could differ slightly from the catalog values used for the the- 
oretical solution. In this melt adapting re-vision an attempt would 
also be made to extend the back focal length slightly to provide more 
working distance between the last surface of the lens barrel and the 
image plane. Even though this represents a difficult parameter to 
move, the performance of the lens is still comfortably within the 
ibO micron blur circle diameter allowance. Therefore, a slight in- 
crease might be possible. 

general conclusions 


To a great extent the success or failure of this lens system will 
be largely dependent upon the reduction of the scattering effects, la- 
tent in the lens system. These effects will be produced by a variety 
of contributing factors as shown below: 

1. Proximity of "ghost" image planes to the true focal 
plane. 

2. Size and quantity of foreign inclusions within the 
glass volume, including bubbles, striae, and Inhomo- 
geneities. 

3. Presence of surface defects such as: 
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a. Sub-surface micro-fractures due to surface grind- 
ing techniques 

b. Surface scratches and digs 

c. Polishing sleeks on the surface 

Non-un i form opt ica 1 coatings including sputters, stains, 
and pin-holes. 

5* Surface dust particles and stains. 

Under the first item above a detailed analysis was performed 
to determine the location of all ghost images due to internally 
reflected energy. Since the first set of energy reflecting pairs is 
most critical It is important that the reflected Image formed by the 
first lens element is not located near the true image surface. As 
shown in the attached computer output the gncst image formed by tiie 
combination 3> 2 is located at 10.021 inches from the last vertex. 

This Is well removed from the image plane and will not contribute 
any spurious image signal to the detector. The next ghost image of 
importance is formed by the pair 4, 3 and is located at 2.112 inches 
from the last vertex. This distance also is well removed from the 
true image surface. All other ghost images are considerably attenu- 
ated due to multiple reflections and would not need to be considered.* 

The remaining factors on the above list pertain directly to the 
optical shop finishing techniqu-s and to the perfection of the optical 
coating processes. Since the methods of optical fabrication technology 

Ref, Discussion on theoretical energy scattering - Dr, A. Quesada 
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represent especially advanced areas of skilled expertise by SORL 
optical shop technicians, the actual optical fabrication should 
present no unusual problem- SORL has achieved considerable suc- 
cess and background knowledge with the production of low scatter 
surfaces for the laser fusion industry using high energy lasers. 
However, the quality control of the coatings will demand painstak- 
ing attention. Inspection equipment must be designed to monitor 
the surfaces after coating and those which do not meet specifica- 
tions imposed must be rejected and re-done. 

The lens system, finally, must be assembled In a true, clean 
room environment where the atmosphere is micro-filtered to remove 
all contaminating dust particles. SORL does maintain an electronic- 
ally filtered clean room for such assembly purposes. 
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S. SCATTERING ANALYSIS 


I Introduction 

Techniques to estimate the fraction of light scattering 

from very rough surfaces have been known for several decades 

(1). They are based on mathematical models that assume the 

phenomenon to be capable of description within the scope of 

geometrical optics, by combining statistical methods with 

estimates of the path differences of- luminous rays that have 
been operated on by the irregularities of the rough surface. 

The experimental evidence that was available showed that, 
although consideration of the magnitude and number of 
irregularities led to a qualitatively acceptable interpre- 
tation of the scattering phenomenon, quantitative predic- 
tions were incorrect, particularly when the spatial distrib- 
ution of scattered energy was desired. The need for closer 
agreement between theory and experiment led the early 
investigators to the conclusion that a theory from which 
acceptable quantitative estimates could be derived required 
that diffraction effects be taken into account. This 
requirement becomes mandatory ^when surface or refractive 
defects become very small and/or are few in number. Attempts 
to develop diffraction oriented theories have not been very 
numerous. To our knowledge, the earliest one is due to A, 
Marechal (2) who treated the problem in the light of the 
theoretical framework of Fourier methods described in detail 
by Duffleux (3)- In the present work we shall extend some 
of the techniques introduced by Marechal to include explicitely 
the effects of off-axis geometry and to take into account 
the effects of microscopic inhomogeneities in refractive 
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elements. In agreement with references (2, and 5), 
we shall be concerned primarily with optical surfaces of good 
quality, whose irregularities are small enough to make all 
deformations imposed upon a wave front small in comparison 
to the wave length. The results to be presented are of 
two types. We shall obtain expressions for total integrated 
energies and, later, we shall develop the expressions for the 
corresponding energy distributions. Detailed derivations 
of equations that are not directly concerned with the 
evaluation of those energy expressions will be relegated to 
appendices. Whenever iL seems appropriate, we have applied 
the equations to the numacical evaluation of scattered 
energy ratios for conditions that may arise in practice, 
or for the establishment of order of magnitude estimates 
that can serve as guidlines to determine whether detailed 
computer evaluation of multiple integrals are necessary 
to improve precision to a satisfactory degree. The appli- 
cability of Fourier methods is a consequence of simpli- 
fying assumptions regarding the geometry of a system. When 
the conditions that allow the simplification are not met, it 
becomes necessary to use a more general formulation in which 
the diffraction integrals have to be evaluated either 
numerically, on high speed computers, or asymptotically, 
by means of methods such as Laplace^s, Kelvin's, or Debye's. 
The particular choice is primarily dependent on the limits 
that appear on the integrals to be evaluated, 
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The next two sections are devoted to the development of the 
theory on the basis of Fourier methods. The extensions to 
the more general formulation is, in most cases, immediate. 
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I I The total energy diffracted 


Suppose we have a refractive surface S, approximately 
flat, that creates small irregular deformations on the 
wave front. To compute the total energy diffracted, we shall 
Imagine that a lens L is placed as in the following figure 
to transform the quasi-plane wave front into a nearly 
spherical front 2^ having the same irregularities. To 
calculate the distribution of amplitudes at a point P (x, y) 
of the focal plane of the lens, we shall use on the wave front 
the angular coordinates ), the angles that the 

direction of the line joining the origin with the point on 
the wave front makes with the planes xO'^ and , 

The use of the angular coordinates simplifies the exponent 
of the phase function that appears in the diffraction 
integrals. Further simplification to transform the 
integrals into Fourier integrals requires that certain 
conditions be satisfied by the coordinates of the obser- 
vation point on the focal plane or near the focal plane 
of L. These are discussed in appendix 2. We know that 
under the conditions that allow writing the integrals as 
Fourier integrals (6) it is possible to express the ampli- 
tude of the electric field E (x, y) at the point P (x, y) 
by Huygens' principle, according to which if F ^ ) is the 
amplitude distribution on a sphere centered at 0, then 
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( 1 ) 

where R is the radius of the reference sphere and k = 

If the wave front is deformed by an amount A with respect 

to the reference sphere the amplitude on the sphere is of the 

r r ‘kA, where E is a constant over the wave front, if 
form E e o 

o 

we now postulate that A Is small in comparison to ^ , 

we can write 



( 2 ) 

We can choose the radius of the reference sphere mh such a 
way that the mean value of A vanishes (A = O) . In these 
circumstances, the amplitude distribution over the aperture 
is the sum of two terms: a) a uniform amplitude 

and b) a fluctuating amplitude of mean value zero ) = 

ik ). 

Consequently, the resultant field E (x, y) is the sum of 
two terms. The field is the field diffracted by the 
aperture utilized. It is very small, except near the 
geometrical image of the source, within a domain whose area 
varies inversely as the area of the aperture. We shall be more 
interested in the field whose amplitude generally 
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decreases far more slowly. The defects represented by the 
function are, in general, small in extent by comparison 
with the extent of the aperture. Consequently, they 
diffract light into an area far more extensive than the 
area over which has appreciable magnitude. To compute 
the energy diffracted, we must introduce a more precise 
definition. We shall adopt, with Marechal , a definition based 
on the integral of the flux due to the amplitude E2. Using 
Parseval's theorem, which expresses energy conservation, we 
can write this energy as proportional to 



( 3 ) 


which, in view of the definition of F is equivalent to 



w 


This integral is recognized as the mean square 
A. It is convenient to introduce the relative 
scattered, by division of the above expression 


2„2 


energy, which equals E^ R 






value of 
energy 
by the total 
. Then , 
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Thus, the energy scattered over the image plane, expressed 
as a fraction of the total energy, is proportional to the mean 
square value of the optical pathvlength error We may 
recall that a similar result is obtained in connection with 
the loss of illumination at the center of a diffraction 
pattern when small aberrations are introduced. It is known 
that in this case, the relative loss of Illumination is 
precisely equal to k^ i and that the energy taken away 
from the central spot is distributed over the focal plane. 
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ill Computation of A for a scratch, a bubble, and micro- 
scop I c i r regu 1 a r 1 1 i es . 

Let us apply Eq, (5) to an irregularity in the shape 
of an ellipse. Limiting cases are clearly a narrow straight 
groove that can simulate a scratch and a crater*^1ike depression 
that may be representative of a flaw in an ant i -ref lection 
coating. In line with the discussion of section 2. 1, the 
dimensions of the ellipse are small in comparison to the 
effective area of the domain of integration. From appendix 
1, we have the A values 

A “ 2ne ws 1 for a reflector, and 


A - (cost - ncosi) e for a refractive 
surface. Since the dimensions are small, the angles i and r 
are essentially constant. Also, strictly speaking, to 
satisfy the condition A = 0, we should modify the defini- 
tions of A above. However, since the area, A, of the i rregu- 
1 ar ity is such that A<<2, where II is the area of the 
aperture, the correction to e is of order A/2 and, there- 


fore, entirely negligible. We then have 


we Lnen nave , ^ \ - 

, rV' /(fTT^vY'' 

d 


f 0 / 
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' Hence, 


P~ (V) A* 

A' 

where is the solid angle subtended by the aperture at the 
focal point. This equation can be modified to read 

p ^ / ^/r \ 

A V X j 

where a. and b are the semi-axes of the elliptical irregularity 
and A is the effective area of the aperture. For a refrac- 
tive Surface, for example, we have in the worst possible case 

A = (1-n) e 

Therefore, 

/X7r\^ 

\t 1 

For an axi-symmetric system whose effective aperture is 
defined by an f-number we can write 


In other words, the fraction of the energy scattered is 
proportional to the area of the irregularity, to the square 

-21- 
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of the depth of the irregularity measured in wave lengths 
and inversely proportional to the effective area of the 
aperture. For a bubble of radius a inbedded in a refrac- 
tive surface, we have approximately. 



and 

P. 






/ ^ 

is 





Finally, for microscopic irregularities, let us consider 

a very simple model in which they are distributed in such 
a way as to cause a random error in the mean position of 

the wave front, for which the probability of finding a 

value of c between e and e & de is given by the expression 





where t represents the thickness of the rej;ractlve element 
and An is a representative value of the index fluctuation 
in the medium. Taking the ensemble average, we find that 
with N the average number of volume elements with distinct 
values of the index encountered by a luminous ray. 
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Hence, 


L (fr 



P ‘■I I’'" 

^ - A/' V A / 

Since An can be as small as 10 , the energy loss from this 
source of irregularities is generally very small, as N is of 
the order of 100 0/cm. 
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1 V Distribution of Intensities. 


To characterize the intensity distribution, ^we must 
construct the function 


D ^s‘i) 



X 




( 6 ) 


The function (x, y) , by virtue of Huygens' principle can 
be expressed as a Fourier transform of the function F2 . 

Since, in general, the function F^ ) is random, it 

may appear that we can only calculate the mean value of its 
square. This, however, is not the case. It Is possible to 
calculate the energy distribution because this energy distrib- 
ution is the Fourier transform of the correlation function 


of F^. Recall that 


E2<x. y) = 'T' 

'' IJ 


if P 


p: (’<,/“) ^ 




f ‘ 


7 


(7) 


2 

The energy distribution is j^E2(x, y) j , a function whose 
Fourier transform may be found by means of the correlation 
function, that is, its own Fourier transform. The deri- 
vation is straight-forward and may be found, for example, 
in reference (7). The result can be stated as follows: 
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\f Q is ) is 

form of g (x, y) , then 


the two“d ) mens i ona 1 


Fourier trans- 



and s imi larly , 




(8b) 


where the script F stand for “the Fourier transform of." 
Let us apply Eq. (8b) to (x, y) , identifying 
and denoting the resultant transform d ( ). We 

have 

-R 

( 9 ) 





Note that in Eq. (9), we have introduced the non-dimen- 
sional displacements ^ is necessary 
to preserve the dimensional homogeneity of the arguments of 
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the functions a (d, g) and A (d -f«d', gfg'). It is 
sometimes convenient to write the correlation integral in 
the form 



( 10 ) 


where is the solid angle subtended by the aperture and 
ip (d‘ , B') is a normalized correlation function 





Using Eq. (10) establishes a direct link between the energy 
distribution, its Fourier transform, and the mean square 
value of the path length errors that we calculated in 
section 2.1. Eq. (10) indicates that if a defect covers 
an appreciable portion of the aperture (for example, an 
error of figure) the correlation function decreases slowly. 

In the case the functions A (d, f>) and A (di^-d', g-f-p*) 
are nearly identical and the integral of the product is 
approximately equal to_SLA? If, on tne other hand, A (d„ g) 
varies rapidly, the two functions contribute to the integral 
only over a limited range of displacements. In either case, 
the energy distribution, | E^ (x, y) | is obtained by inversion 
of the transform, Eq. (10). That is, 
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or equivalently, ^ y ^ ^ - Z??'. J J,^j 


D ov:i .-- 




J j ‘ i ^ 


r- H'l; ''!< ''t,^n cJ (4 , H- 

<// X ^ .A -'X y 


A 


( 12 ) 


•J ( f , j( 


wliere ( ' <^ ^ K ) denotes the double integral Inverse 

A * ■/ 

transform. 

This expression suggests a simple way of measuring the 

energy distribution. Suppose we can determine the mean 

square value of A and can measure the ratio of the flux 

2 2 

D (x, y) to the total flux (-R ) at given locations, 


but as functions of the wave length. Suppose that 


subsequently we construct the curve representing 
£/‘(x, y) as a function of x/^ and y/A. According to Eq. (12), 
this curve represents the product A p. Since A is a very 
weak function of the wave length, upon division by M 
and subjecting the function ^ to the condition that its 
Fourier transform is equal to unity when(?( = P=0, there 
suits the function ^ j . The spatial energy 


resi 


distribution for a given wave length follows immediately 

by a simple change of variable. To avoid interference from 

the intense direct images, the experiment should be conducted 
with a black hole to eliminate the contribution from the 

direct images. The construction of such a device is des- 
cribed in reference (8). 
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The decrease in illumination is again a Gaussian function^ 
Again, these examples show that large scale defects, such as 
figure errors, lead to large values of oj and, consequently, 
energy distributions that fall off rapidly. On the contrary, 
when the defects are small, the correlation has short range 
and the diffracted light is spread out widely over the 
focal plane. 
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V Examples of Energy Distribution . 

Let us consider the case of an irregularity of rectan- 
gular shape, of dimensions w and A, simulating a small 
scratch. In this simple case, the normalized correlation 
function Is 





and, therefore, 



As a second example, let us consider an autocorrelation 
function of Gaussian character. VJe might expect such a 
correlation function to describe the wave front distoi — 
tions introduced by microscopic inhomogeneities. Let, 
then , 




f 




to 


where w is a parameter that fixes the scale of correlation 
and, in general, can be interpreted as indicating that the 
irregularities of the wave front are seen from the origin 0 
as subtending an angle of order w. The function 
is in this case ^ ^ V ' 

'A.' / 
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VI SUMMARY OF SCATTERING ANALYSIS 


The theory that was developed to <jeal with levels and 
spatial distribution of scattered light can be applied with 
little modification to define limits on the number and extent 
of scattering centers that may be tolerated in the tracker optics. 

It is important to bear in mind that a theoretical study can very 
effectively provide tolerance criteria to be met by the optica] 
designer, rather than predict in detail how a system will perform. 

The latter requires that a system already in existence be examined 
to determine the type, magnitude and distribution of bulk and sur- 
face defects in each component. if this is known then the theory 
can predict how much light will be scattered and how will the energy 
be distributed. Theoretical considerations are also useful in in- 
terpreting the results of scatter measurements on a system being 
tested. This is particularly true of the stretch correlation, that 
provides outputs that are invariant under changes of scale. 

When the theory is applied to determine tolerances it becomes 
apparent that the most critical parts of the system are the front 
elements. This is because the intensity of the scattered light reach- 
ing the image plane is proportional to the flux incident on the irreg- 
ularities responsible for diffracting the light and introducing dis- 
tortion of the wavefront. The effect of the irregularity itself, as 
far as diffraction effects are concerned, is the same regardless of 
where that irregularity is located. It leads to distortions of the 
wavefront characteristic of the type and dimensions of the irregularity. 
A small bubble, for example, that diffracts 1/10 of ]% of the light 
may be tolerated in an inner element shielded by geometry from direct 

sun light at all angles of interest. The same bubble would totally 
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destroy the system performance if present in the front element 
under conditions that allow exposure to direct sun light. In 
this case, even the fact that the distribution of scattered 
light will in general be oriented at some angle with respect 
to the axis of the system would not be sufficient to compensate 
for the direct effect of incident flux. The incident flux is gen- 
erally strong enough to cause an intolerable loss of contrast in 

spite of the small fraction of side scattered light that reaches 
the image plane. 

When the equations were specialized to account for effects of 
small bubbles it was found that the front elements have to be ess- 
entially bubble free. Although the equations were used to define 
some limiting upper size for a bubble, they may be interpreted also 
as indicators of how many smaller bubbles can be allowed. The trade- 
off is not immediate because it must be recalled that even though a 
small bubble is responsible for a small distortion of the wavefront, 
the light that it scatters in the forward direction is spread over a 
much wider solid angle that would be the case if the size of the 
bubble is allowed to increase. Reference to the table on page ]lf 
is necessary to estimate how the angular factors are to be changed 
when the bubble size is decreased. 

A similar argument can be used to estimate the effects of sur- 
face defects. If the equations applicable to a scratch are used, we 
find that in the worst case no scratch of depth of order of the wave- 
length and longer than a few millimeters can be tolerated for the front 
surface or any surface that may be illuminated directly by the sun. 
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The worst case corresponds to a scratch oriented in a direction 
that will cause maximum spread of the diffracted light. It is im- 
portant to remember that the quality of the surface finish or the 
quality of a coating can be evaluated in a variety of ways. For 
the small defects that we have considered electron microscopy, 
phase microscopy, dark field light microscopy and ellipsometry 
are suitable. For larger surface areas phase contrast and knife 
edge tests can be used. The measurements by any of these techni- 
ques can be used as inputs to the formulation that has been devel- 
oped, to confirm that irregularities have been kept below the limits 
needed to meet performance spec i f i cat ions. 

From the literature references, these are A main points which 
require special attention in terms of securing low scattering in 
lens systems. They are: 

1. For a multiple lens element system the effect of 
flare is less at higher relative apertures than 
lower apertures. 

2. For extended self-luminous objects, such as the 
horizon, the magnitude of the flare is propor- 
tional to the number of air-glass interfaces and 
not determined by the relative aperture. 

3 . In designing lens systems for low scatter applica- 
tions, ghost reflected images adjacent to the Image 
plane must be scrupulously avoided. 

A, Good optical coatings do reduce the absolute magni- 
tude of scattered energy, but porous films will 

slightly increase the amount of scattered energy. 
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VI I EFFECT OF BUBBLE IN FRONT ELEMENT, DIRECTLY ILLUMINATED BY SUN 


Eg^s. (3) and (6) of Part I can be used to determine the 
effect of a small bubble in the front element of the proposed 
systems, when illuminated by direct sunlight falling on that 
element at an angle beyond the field of view. Several simpli- 
fying assumptions will be introduced, but in every case the sim- 
plification will be chosen in such a way as to produce 'worst* 
effects. This way we obtain criteria that, if met, will insure 
proper performance of the system. 

Since the correlation function will in every case be analo- 
gous to the first example in section 3-2, the worst case can be 
approximated by distributing the total energy evaluated from eg. 

(3) over solid angles bounded by annular areas of radii equal to 
the distances corresponding to consecutive zeros of the function 

sin X R . We shall identify co with the radius of 

the bubble. It is convenient to express the argument of this func- 
tion in terms of angular quantities measured from the direction of 
the incident light. The zeros, then, occur at angles 0^ such that. 





TABLE III. 


co/X 


^2 

®3 

'4 

0.1 

84.3° 

87 . 1 ° 

88.1° 

88.6° 

0.5 ' 

63.4 

76 

80.5 

82.9 

1 .0 

45 

63.4 

71 .6 

76 

1.5 

33.7 

53.1 

63.4 

69 

2.0 

26.6 

45 

56.3 

63.4 

3.0 

18.4 

33.7 

45 

53.1 

5.0 

11.3 

21 .8 

31 

38.7 

10.0 

5.71 

11.3 

16.7 

21.8 

20.0 

2.9 

5.71 

8.5 

11.3 

50.0 

1.15 

2.29 

3.43 

5 . 7 ] 

100.0 

0.57 

1.15 

1 .71 

2.29 

200.0 

0.29 

0.57 

.86 

1.15 


The normalized peaks of energy density within each one of the areas 
bounded by zeros are 1, 0.0^7* 0 . 003 ^ and 0.002A, respectively. 


From eg. (3), we find that the total energy scattered by the 
bubble is, 

E„ = C 

D O 

2 

where E is the amplitude of the E-field and CE can be taken as the 
o o 

energy density of the incident light. 

Now, according to eg. (3) and the table shown before, the smaller 
bubbles scatter little light, but the scattering Is over a very wide 
angle. On the other hand, the large bubbles scatter far more energy, 
—2 

(large A ) but that energy is concentrated in the forward direction 
and relatively little Is deflected to the sides. For example, when 
the bubble size is of the order of lOX, the first zone extends to 
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less than 6*^ from the direction of incidence. Letting f be the 
attenuation factor due to angular effects we can then write: 


P = 


Energy density of scattered light in direction of g-th maximum 
Energy density of incident light 


2 7T 


A ^ 


(2) 


If p is to be less than some small fraction E we must have 


1 Stt 


6 2 Tlo 

oj A 2 


f < E 


1/6 


or 


16 TT f 


R 

o 

A (n-1) 


1/3 


(3) 


The light that Is scattered is allowed to pass through other 
elements of the system. When the angle of scattering is such that the 
fraction f reaches the focal plane, a further weakening of intensity 
takes place because the front element in general will not be In focus. 

The decrease in Intensity can be estimated when the position of 

the image plane for the bubble is known. The problem is similar to 

that encountered when one computes the flare of a system due to second-* 

ary reflections of air-glass interfaces. In this case, it Is known* 

that the amount of flare light is proportional to the square of the f/ 

2 

number of the system and is inversely proportional to the (distance) 
See, for example, Kuwabara, G., J. 0. S. A,, 43, 53 (1953)- 
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of the secondary image from the focal plane. The latter effect 
is more serious and in designing an optical system it is desir~ 
able to avoid any pair of reflecting surfaces which produce a 
reflection image in the vicinity of the image plane. For the 
6 -element system, based on the figures obtained for the position 
of secondary image planes at normal incidence, we estimate that 
the reduction is of the order of 10. In general, we can call 
this factor f^ and write instead of eg. ( 3 ) 



~k c 

If X = .5 X 10 cm, = K5, = 3 cm, E = 10 “^, 

f = 0.001 , fj =0.1 

0) < 'v 20 

A 

This means that bubbles in the first 2 elements must be essen- 
tially absent. The last condition may be relaxed somewhat, by re- 
calling that the derivation of eg's. ( 3 ) and ( 6 ) required the wave- 
front distortion to be small. For bubbles 'vlO A this is not strictly 
true and a more geometrical approach indicates that most of the scat- 
tering is forward, with edge effects superimposed on it, and amounting 
to a fraction proportional to some effective volume of the edge as com- 
pared with the entire bubble. 

Equally important is the problem associated with ant i -ref lect ion 
coatings for the first 2 elements. Coatings are effective in reducing 
the flare caused by reflected light, but one must be very careful to 
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have coatings very homogeneous. Porous films enhance slightly 
the scattered light and tend to reduce contrast in the vicinity 
of a bright image. 
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VIII COMPARISON OF 6-ELE/£NT AND A-ELEMENT DESIGNS 


Since a fractional rejection level of 10 ^ will 
attennate sunlight to the equivalent of a uniform back- 
ground brightness of 0.1 foot- 1 ambert , we shall examine 
the expected performance of the systems in terms of fractional 
rejection levels. The same criterion will be used for 
other sources such as the full moon or the earth*s sun- 
lit horizon. For conversion, it is convenient to use the 
fact that a star of visual magnitude 0, corresponds to 

illumination outside the earth's atmosphere of 2.65 x 10 
2 

lumens/cm . Thus, the full moon, with an equivalent 
visual magnitude of -12-70, leads to an illumination of 
3.2 X 10 ^ lumens/cm^. The effect of the full moon, as 

far as incident flux on irregularities is concerned, is 
then nearly 300,000 weaker.** 

The computations that follow were carried out under the 
assumption that glasses will be specified, at least for 
the first two elements of each system, which will have: 

a) Homogeneity: H2, that Is, An = ± 5 x 1 0 ^ 

b) Bubble Group 0, that is, total bubble cross- 

2 

section in mm /lOOcc volume of glass be- 
tween 0 and 0.029. For the calculation, 
we require selection of the glass to have no 
bubble >*60>j and it will be assumed that 
smaller bubbles that may be present will have 


**Allen, C. W,, Astrophys i cal Quantities, U. of London, 
Athlone Press (I963) (pg. 1^5)- 
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2 

a maximum cross section of 0.015 mm /lOOcc with a mean 
size approximately 20X. 

The volumes of the first two elements for the two 
systems were estimated to be 30 cc and 27 cc for the 6- 
element system and 9.1 and y.Zcc for the ^-element system. 

For these circumstances, we can apply the equations 
developed in Part I of the theoretical analysis and 
construct Table I for the cases where the sources are 
outside the field of view. 

These levels are comparable to the total values measured 

for good coronograph lenses where (for those lenses not limited 

by florescence) the fraction of solar light scattering was 
-6 -6 

between 1.2 x 10 and 5.7 x 10 (See Newkirk and Bohlin, 

App. Optics, 2, 131 , 1963 ). 

For the case of the full moon within the field of 
view, we have similar results. The most serious effect 
is due to the small angular deviations that we must now 
consider and the possibility of involving all elements of 
tile system. 

Since, however, the inci'dent flux due to the full 

moon is some 300,000 weaker, in spite of the requirement 

for the equivalent background brightness to be 0.05 

foot-1 ambert , the changes in the factors shown in Table II 

5 

as compared with Table I are well v/ithin the 10 factor 
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^ that the weaker source and equivalent background allow. 

Thus, either system will satisfy this specification. 

Again, we emphasize that the figures represent worst cases. 

A similar analysis for the case of the earth‘s sunlit 
horizon within the field of view shows that the specifica- 
tion can be met at angles of about 3— away from the edge of 
the horizon image. At closer angles, the presence of 
bubbles in the size range 20X to 30A may lead to areas 
where the specification is violated. 

Ray tracing was used to determine the location of 
secondary image planes and to estimate the effects of the 
secondary images as sources of unwanted stray light. 

Assuming the reflectivity of air-glass interfaces to be of 
the order of 1^, we can summarize the results in Table Ml. 

When the full moon is within the field of view spurious, 
images may be formed with illuminations in the image plane 
as bright as 6.4 c 10 ^ lumens/cm^ and 1.6 c 10 ^ lumens/cnri^ 

(for 6-element and 4-element systems, respectively). 

A fourth magnitude star having its energy distrib- 
uted within 12arc minutes square field of view corresponds 

— 6 2 

to an i 11 uni i nat i on of 2 x 10 iumens/cm for the 6-element 

_7 2 

system and .8 x 10 lumens/cm for the A-element system. 

-ho- 
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‘We, therefore, have signal to noise ration of 


For the 6-element system: 


S/N 


2 X 10“^ 

6. A X 10 


30 


For the 4-element system: 


S/N = 


.8 X i0“^ 

1.6 X lO'^ 


1/2 


which shows the 6-element system as a more promising 
system in this respect also. 
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TABLE I 


MAXIMUM FRACTION OF INCIDENT INTENSITY SCATTERED 


SYSTEM 1 

ANGLE 

DUE TO 
BUBBLES 

DUE TO BULK ( 

INHOMOGENEITIES 1 

DUE TO SURFACE 
DEFECTS* (PER DEFECT) 




i 


6-E 1 emen t 

30- 

1.6 X lO"^ 

3 X 

,0- 1 


1 

o 

X 



7.5 X 10“^ 

lo"^ 

4-E 1 ement 

30- 

1 

is X lo'^ 

1.5 X lO"^ 

lo'^ 

1 

1 

16^ 

1.2 X lO"^ 

1 . 

k X 10~^ 

|lu'^ 


A EACH DEFECT EQUIVALENT TO A MiCROSCRATCH OF DEPTH A, 
WIDTH == 1 MICRON AND LENGTH = 1000 MICRONS IN WORST 
POSSIBLE ORIENTATION (i.e. ANGULAR FACTOR = 1). 


_/, 3 _ 



TABLE I I 


MAXIMUM FRACTION OF INCIDENT INTENSITY SCATTERED 


SYSTEM 

DUE TO 1 

BUBBLES 1 

DUE TO BULK 
INHOMOGENEITIES 

DUE TO SURFACE 
DEFECTS (PER DEFECT) 

1 

6“*E 1 ement 

1 

1 

X 

o 

1 

o 

X 

lo"^ 

|4-E] ement ! 

1 j 

io“^ 

1 io“^ 

i 

j 

10-6 

i 


TABLE I I I 

SPURIOUS IMAGE REJECTION 


SYSTEM 

AT lb- 

AT 30- 

AXIALLY 





1' 

■6-El ement 

Rays do not 
reach image 
plane 

Rays do not 
reach image 
plane 

< Z X lO"^ 

for reflections due 

to 1st two elements 

^"E 1 ement 

< lo"^ 

1 

O 

V 

< 5 X lo‘^ 

for reflections due 

to 1st two elements 
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C. Fabrication Details 


I . General Fabrication Summary 

The complete specLfl cat ions covering the optical 
fabrication tolerances are itemized in Table I* Most 
of these values are rather modest and easily held. 

However, some of the tolerances pertari.ning to 
individual element wedge, and leas centering error, 
are severe and must be carefully controlled. There is 
also one air space tolerance value which may require 
special assembly adjustment in order to obtain peak 
theoret i cal performance. 

I I Glass Material Selection 

The particular glass materials chosen for this design were 

dictated by tiie design requirements discussed in the 

previous sections. However, the quality is determined 

by specifying “Schl ieren Grade‘S, PH/jwith Special 

Bubble Selection. This insures that the variation in 

“6 

i nhemogenei ty will be of the order of ^ 1 x 10 with- 
in each blank, stress birefringence ^3"'^ nm/cm and 
maximum inclusions (i.e. bubbles) of 0 - ,023 cross 
section per 100 cc. 

I 1 ( Distortion Characteristics 

The general distortion characteristics of the lens 
system are summarized in Table II. These values were 
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obtained directly from the computer print-out of spot 
diagrams giving the image centroid heights. The tolerance 
V alue of 3arc minutes for distortion was then converted 
directly into the allowed linear displacement value 
of O.OA^mm. Then, the calculation of the residual 
distortion in the lens system was derived directly from 
the value in the table. 

I V Transmission Characteristics 

Table III shows the given catalog values of nominal 
Internal transmittance values for the selected glass 
types. These glasses all show high values of trans- 
mittance at 400 nm, which is primarily why they were 
chosen. A sampling of various firms specializing in 
optical coatings also indicate that the nominal value 
of 0 . 5 ^ represents the lowest- ref lectivi ty value to be 
expected over the range from 0.40 to 0.62 microns. 

These coatings will satisfy MiL-C -675 specifications for 
abrasive resistance and weathering characteristics. 

V Mechanical Considerations 

The cell will be made of a stainless steel which most 
nearly matches the glasses with regards to the thermal 
co-efficient of expansion. The elements will be edge 
• shimmed at assembly for zero clearance. A window, for 
protection of the optical elements, is provided at the 
front of the cell and is removable by means of one 
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retaining ring. 

The lens will incorporate a lockable focussing arrange- 
ment providing 2mm of motion on each side of best focus. 
The focus housing will engage the 2.118-48 thread and 
a shoulder will bottom on the system face. The lens 
cell is piloted in the focus housing on two lands and 
moves axially by means of a 12 pitch thread. Locking is 
accomplished by a set screw. 

Light baffling will be accomplished by providing all 
retainers and spacers with deep grooves such that light 
entering the grooves cannot continue downstream except 
after two or more reflections. Ail retainers and spacers 
will be surface treated for minimum reflectivity with 
minimum specular component. 
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TABLE I 


TABLE OF TOLERANCES 


SURFACE AR 

SURFACE 

FIGURE 

U 2 

±0.014 

.50A 

3 

±0.012 

.50A 

2{4 

±0.003 

• 33A 

5 

±0.015 

• 33A 

3^6 

O 

o 

o 

+1 

. 25 !^ 

7 

±0.001 

.25A 

-8 

- 


4(9 

±.012 

.25A 

IG 

± . 002 

.25A 

5{11 

±.002 

.25A 

12 

±.001 

• 33A 

8{13 

±.005 

• 50A 

14 

±.246 

.50A 

Note : 

All glasses must be Schlieren 


-Special Bubble 

Selector. 

*MAY 

REQUIRE ASSEMBLY 

ADJUSTMENT 


(ELEMENT CT) 
Ad 

AIR SPACE 
As 

ana 

±.0040 

±0.0014 

±0.0016 

±.0011 

±0.0018 

± 0.0008 

±.0015 

±0.0012 

±0.001 1 


±0.0012 


±.0014 

±0.0006* 

±.0020 

±.00)7 

±0.0031 

+.0006 

±.004 


±.0016 

Grade, PH^ quality witii 



General Surface Frni: 


Element I and I I 
Elements III, IV, V 
VI 


WEDGE 

ELEMENT TILT 
LENS SEAT 

ELEMENT 

DECENTER 

0.0012" 

0.0011" 

0.0010" 

0.0004" 

0.0004" 

0.0008" 

0.0002" 

0.0003" 

0.0004" 

0.0004" 

0.0002" 

0.0006" 

0.0002" 

0.0003" 

0.0004" 

0.0005" 

0 . 0007 " 

o.UOly" 


Target 

Requ i red 

0/0 

5/0 

5/5 

10/5 
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table I i 



DISTORTION 

VALUES 


FIELD ANGLE 

GAUSSIAN HEIGHT 

IMAGE CENTROID 

DISTORTION* 

(Rad i us) 

(mm) 

(mm) 

(min) 

0.000 

0.0000 

0.000 

0.00 

0.026105 

1.3057 

1.2992 

0.44 

0.0522! 

2.6141 

2.6043 

0.67 

0.078315 

3.9278 

3.9128 

1.02 

0.104420 

5.2497 

5.2268 

J.56 

0.130525 

6.5826 

6.5558 

1 .83 


Tolerance Value = 3-0 min. 
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TABLE I I 1 


TRANSMISSION & REFLECTION LOSSES: kOO NH 


MATERIAL 

INTERNAL 

TRANSMITTANCE 


5mm 

25mm 

LAK-8 

0.988 

0.940 

SF-19 

0.985 

0.920 

SF-1 

0.980 

0.890 


Reflection losses for 12 air-glass interfaces, assuming a very 
realistic coating value loss of 0.5%, will amount to 5.8!^. 

The total absorbtion losses will be, nominally, less than 
15.1^. Therefore, the average transmittance of the system 
will be about 79.1^. This value is somewhat pessimistic 
since the absorbtion values are given for plane parallel plates 
of glass and the lens elements are much thinner at the edges. 

The coatings will be broad band coatings, peaked for 4,0008. 

The target value of ^0% than appears to be definitely a practical 
working value. 
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APPENDICES 


Appendix 1 . Calculation of Optical Path Errors Due to 
Surface Defects. 


Consider first a refracting surface separating two 
media of indices n^ and n^. If a defect of depth e is 
present, the difference in optical path for the ray that is 
refracted at I and the ray that would have been refracted 
at I' if the defect were absent is, from the figure, 


X 











in 






p 

p 


Al-AI ^T/\ 




/ 


Let 


TP-Tf-' = x£‘- rip , 
L = t), (aI'-az') -f H, (rP-T'^ (••') 


4 Iw = vi, \ ,* n. I'f-i: ' - ) 






then. 


^ 


II ii: -N 
Pi P “• jyp'- ^ 


(a. 2) 


where N is the unit normal to the surface in the absence of 
defects. Substituting Eqs. (a. 2) into Eq. (a.l) leads to 
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(a. 3) 


Since e is assumed to be very small, 3 ~ 6‘ and 


i: P - -TP ~ - L'/Jf -■ pj‘] 


Hence , 


)'Il 


The law of refraction in vector form reads, 


- nj i = (n 2 cos r_ -/J^cos i ) N 


(a. 4) 


Therefore, we can write 




But R • The deformation of the wave front is 

A = -^L. 

Consequently, 


A = (n^ cos i - n^ 
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For a reflecting system, we have 

A 

F ' p 

S'L, - iTi XT ' ^ ] 

But the law of reflection in vector form reads 

^ y\ 

r = - i +2 cos i N 

Hence , 

to. sf'^J 

Again, N Hence, the deformation of the wave front 

i s 

A = 2 n e cos i 

(a. 6) 
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_Appendix 2. Computation of the Distance That Appear s in the 
Diffraction Integral. ^ 


If a converging wave front has phases represented by 
a phase function F (x‘ , y', z'), the amp i I tude of the disturb- 
ance that reaches a point P (x,y,z) is given by the classical 
Huygens- Fresne 1 formula (see, for example, ref. 6). 





(b.l) 


To evaluate d^, we consider an axially symmetric (or nearly 
so) wave front whose vertex is situated at a distance f from 
the origin of a coordinate system whose z-axis is in the 
direction of the symmetry axis. Let R represent the 
radius of curvature of the wave front at the vertex. From the 
figure it is clear that the distance from a point M on the 
^ wavefront to an 


A 


p(xV/.i) 


r 



\ 


-5J,- 
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observation point near the plane x-y is 





(b.2) 


The position of the point M can be more conveniently 
specified by the coordinates r and (j> plus the distance 
along the axis e , 

V/e can then write for 



-f 


^.i') - ij ^ 


(b.3) 


For the special but important cases where the wave front 


is spherical or 

paraboloidal, we have that 


H i,'- 


y L 

and 

■ X (\ 


respectively. 

A simplification is obtained when all five of the conditions 


hold y- , -•! 

V.' - J. 

i- ^ 




T 






J. 
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In this case 



This is the usual form of for diffraction work. We 
shall be interested in a different form and aiso in the 
case where z — 0, that is, we want to have the observation 
point in a plane, that in general will be the focal 
plane of a lens or a mirror. Let us introduce the new 


variables . 

"4? ' 




As seen from the following figure, the angl es^and 
are in orthogonal planes and allow us to write d 2 in 
the form 


(b.5) 

Since the condition r/f « 1 implies that the angles 
cL and ^ are small, we can replace the sines by the 
angles and have only the linear combination as the 
variable part of the exponent, in which case the integral 
becomes a Fourier integral, with exponential term 
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Appendix 3 « Sample l»cattering Computation . 

Computations are based on the following three equations that 
represent the amount of light scattered by an irregularity 
as a fraction of the flux incident on it. 


For Bubbles: 


e = 1677 { A (n-1) f (8) f (d) 

Ro 12 

For scratclies ; 

e * (27 t) 3 (n-1)^ (h)^ Ax fiy f (e) f (d) 

J 


For bulk i nhomogene i t i es : 
A N ^ 


Ro 


(d) 


where 


A = wave length of light 
n = index of refraction of glass 

Ro “ radius of curvature of wave front forming image 
0 ) = radius of bubble 

Ax, Ay = dimension of scratch determining projected area 
h - depth of scratch 

t — thickness of lens or optical element 
An = fluctuation of index in optical element 
N — av. number of doma i ns/un i t thickness showing the 
fluctuation An* 

” factor to account for angylar distribu- 
tion of scattered 1 ight 

factor to account for defocussing effects 
that tend to spread out the light over the 
image pane. 
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• For bubbles, for example, we have selected a glass 

belonging to bubble group 0, specified as having a total 

bubble cross section between 0 and 0.029 mm /lOOcc volume of 

glass. We shall assume a value of 0.010, for 1st element. 

For the 6-element system, we further require that there be 

no bubble or inclusion of diameter 90u or greater. This is 

consistent with the lower limit of bubble group 0. Since the 

1st element of the system has an approximate volume of 30 cc, 

we have a total area, according to specifications, between 0 
2 

and 0.01 (mm) . Assuming a log-normal distribution for the 
bubble present, the mean bubble radius is in the vicinity 
of 20A. The radius of curvature of the front surface will be 
taken as 10cm and for a 20^ bubble with direct illumination 
from a direction making an angle of 30 ^ with the optic axis 

-3 

we have an angular factor *'10 . With these values, and recalling 

that the true image plane is about 2.5 cm behind the last 
element we have 


e = 16 it {.16 


io“^} 

¥ X To o 



r3 


= 1 6it 


(16 X 10 


-10 


) (10 


12 


rs— 

A X 


) X lo'^ 

.3 


= 10 


-4 


1.6 X 10 


-6 
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For this case, we have taken f^(d)==l, to consider 
what might be regarded as a worse case. For scratches, we 
have taken the angular factor to be 1 , since there may be 
cases where the orientation of a very narrow scratch spreads 
the light over a considerable angle from the direction of 
the incident 1 ight. 

For bulk inhomogeneities, we consider an element 
l4mm thick made of galss in homogeneity group H i.An=±5xlO }. 
Since N is of order (lO-^) 'em we have 

e = (5 X ( 27t X ] A i0~^ 

k X 10“^ 1.4 X 10^ X 2.5 

= 1/3 X 10“^ 

also have assumed that when illumination is such that 
more than one element receives light from the source, the 
net effect is additive. 
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.Appendix h . Computer Design Data. 
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